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Three events for the decay Jf"*" — + ■k'^vV have been observed in the pion momentum region 
below the K'^ — > 7r"'"7r° peak, 140 < Ptt < 199 MeV/c, with an estimated background of 
0.93 ± 0.17(stat.)^Q |4(syst.) events. Combining this observation with previously reported results 
yields a branching ratio of -r^vu) — (1. 73^105) x 10"^" consistent with the standard model 

prediction. 



PACS numbers: 13.20.-v, 12.15.Hh 

The decay if + —>■ ti^vv is among a handful of hadronic 
processes for which the decay rate can be accurately pre- 
dicted in the standard model (SM) owing to knowledge of 
the transition matrix element from similar processes and 
minimal long-distance effects [H |3| • The small predicted 
branching ratio, B{K+ tt+vv) = (0.85 ± 0.07) x 10~i° 
3, and the fact that this decay is a flavor-changing 
neutral current process makes it a sensitive probe of a 
wide range of new physics effects Previous stud- 
ies of this decay by experiment E787 at Brookhaven 
National Laboratory and its upgraded extension E949 
have measured B{K+ tt+vv) = (1.47+Jjg) x 10"!° 
based on the observation of three events in a sample 
of 7.7 x 10"'^^ decays at rest with a total expected 
background of 0.44 ± 0.05 events in the pion momentum 
region 211 < < 229 MeV/c above the K+ tt+tt" 
{K^2) peak (pnnl) E787 set a consistent limit of 

< 22 X 10~^" at 90% C.L. based on one candidate in a 



sample of 1.7 x 10^^ stopped if + decays with an expected 
background of 1.22±0.24 events in the momentum region 
140 < < 195 MeV/c below the K^2 peak (pnn2) [alzl- 

In this Letter we report the results of a search for 
isr+ TT^vv below the if^2 peak (pnn2) using 1.7 x 10"^^ 
stopped K'^ decays obtained with E949 as well as the fi- 
nal results on B{K^ ■n'^vv) from E949 data combined 
with E787 data. 

Identification of if + n^i/D decays relies on detection 
of an incoming kaon, its decay at rest and an outgoing 
pion with no coincident detector activity. The E949 ap- 
paratus and analysis of the data in the pnnl region have 
been described elsewhere In this letter, we empha- 
size the apparatus and analysis features most relevant for 
pnn2. 

Incoming kaons were identified by a Cerenkov counter 
and two proportional wire chambers before being slowed 
by an 11.1 cm thick BeO degrader and an active degrader. 
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passing through a beam hodoscope and stopping in the 
scintillating fiber target. Typically 1.6 x 10® K'^ /s en- 
tered the target during a 2.2 s spill with a /tt^ ratio 
of 3. The active degrader had 39 copper disks (2.2 mm 
thick) interleaved with 40 layers of 2 mm plastic scintil- 
lator divided into 12 azimuthal segments. Scintillation 
light from each segment was transported via wavelength 
shifting fibers to a photomultiplier tube (PMT) that was 
read out by time-to-digital convertors (TDCs), analog- 
to-digital convertors (ADCs) and GaAs CCD waveform 
digitizers (CCDs) sampling at 500 MHz [8j]. The active 
degrader was capable of providing measurements of the 
incoming beam particle and activity concident with 
decay in the target. The target consisted of 413 scintil- 
lating fibers (5 mm square and 3.1 m long) packed into a 
12 cm diameter cylinder. Each 5 mm fiber was connected 
to a PMT and read out by TDCs, ADCs and CCDs in or- 
der to record activity in the target coincident with both 
the incoming kaon and the outgoing pion. 

The momentum and trajectory of the outgoing tt^ were 
measured in a drift chamber . The outgoing pion came 
to rest in a range stack of 19 layers of plastic scintillator 
with 24 segments in azimuth. PMTs on each end of the 
scintillator were read out by TDCs, ADCs and 500-MHz 
transient digitizers [lol | and enabled measurement of the 
pion range (i?^) and kinetic energy (E't) as well as the 
7r+ ^ /i+ ^ e+ decay sequence. 

The barrel veto calorimeters of 16.6 radiation lengths 
(r.l.) at normal incidence provided photon detection over 
2/3 of 47r sr solid angle. Photon detection over the re- 
maining 1 /3 of 47r sr solid angle was provided by a variety 
of calorimeters in the region from 10° to 45° of the beam 
axis with a total thickness from 7 to 15 r.l. [HI. [Ill [l^. . 
More extensive use was made by this analysis than the 
pnnl analysis of the photon detection capabilities of the 
active degrader (6.1 r.l.) and the target (7.3 r.l.) that 
occupied the region within 10° of the beam axis. 

This pnn2 analysis was able to increase the signal ac- 
ceptance by 40% and maintain the same background rate 
per stopped as the previous analysis thanks to im- 
proved background rejection primarily due to the addi- 
tion of the active degrader and augmentation of the bar- 
rel veto by 2.3 r.l. for E949. In addition, the improved 
knowledge of the background contributions allowed the 
signal region to be divided into nine sub-regions ( "cells" ), 
with relative signal-to-background levels differing by a 



factor of 4, that were used in the likelihood method 
to determine B{K^ tt^vv). 

To avoid a possible bias, we employed a "blind anal- 
ysis" technique Q in which the signal region was not 
examined until all selection criteria ("cuts") for signal 
had been established, the estimates of all backgrounds 
completed and acceptance of all cells determined. Two 
uncorrelated cuts with significant rejection were devel- 
oped for most backgrounds. After imposing basic event 
quality cuts, inversion of one of the pair of cuts could then 



Process 


Background events 


KTr2 target-scatter 
KTr2 range-stack-scatter 

Charge-exchange 

Muon 

Beam 


0.619 ±0. 150^0 10; 
0.030 ± 0.005 ± 0.004 
0.076 ±0.007 ±0.006 
0.176 ± 0.072lo.?24 
0.013 ±0.013+0-003 
0.011 ±0.011 
0.001 ± 0.001 


Total 


0.927 ±0.168i;j-^3i; 



TABLE I: Summary of the estimated number of events in 
the signal region from each background component. Each 
component is described in the text. 



be used to select a background-enriched data sample con- 
taining N events. Inversion of the complementary cut se- 
lected a data sample on which the rejection TZ of the first 
cut could be measured. The background was estimated 
as N/{TZ — 1). We ensured unbiased background esti- 
mates by dividing the data into one-third and two-thirds 
samples chosen uniformly from the entire data set. Selec- 
tion criteria were determined with the one-third sample 
and background levels were measured from the two-thirds 
sample. In contrast to the analysis of the pnnl region, 
some backgrounds did not have sufficiently distinct char- 
acteristics to permit isolation by cut inversion of a pure 
background sample and permit a measurement of TZ with 
the data. For these backgrounds, TZ was estimated with 
simulated data as described below. 

Table [J summarizes the estimated background levels. 
The largest background was due to /C7r2 decays in which 
the 7r+ scatters in the target, losing energy and obscur- 
ing the directional correlation with the photons from the 
7r° decay that would otherwise be detected in the barrel 
veto. Two cuts that suppressed this background were 
1) identification of 7r+ scattering and 2) detection of the 
photons from 7r° decay. Pion scattering was identified 
by kinks in the pattern of target fibers attributed to the 
pion, by tracks that did not point back to the fiber con- 
taining the decay, by energy deposits inconsistent 
with an outgoing pion or by unexpected energy deposits 
at the time of the pion in fibers traversed by the kaon. 
The target pulse-shape cut identified the latter signature 
by performing a least-squares fit to the CCD samples to 
identify the pulses due to activity coincident with the 
kaon or pion Q. The uncertainty in the Kt^2 target- 
scatter background had comparable statistical and sys- 
tematic contributions. The systematic uncertainty was 
determined by the range of photon veto rejection values 
measured on samples of Kt^2 scatter events selected by 
different scattering signatures in the target or in differ- 
ent 7r+ kinematic regions TH]. There was also a much 
smaller background from KTr2 due to scattering in the 
range stack that was similarly identified by the energy 
deposits and pattern of range stack counters attributed 
to the track. 
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Additional backgrounds included K'^ 



and — > 7r°//+i^ (muon), scattered beam pious (beam) 
and 7r+ resulting from charge-exchange reactions 
dominated by n^fi^u. Simulated data were used 

to estimate the rejection TZ of the cuts that suppress Ke4, 
Kt^2-i and charge-exchange backgrounds. The K^^i and 
Kt^2-i backgrounds could not be distinguished from the 
larger iir^2-scatter background based solely on the 7r+ 
track, and it was not possible to isolate a sufficiently 
pure, statistically significant sample of charge- exchange 
events on which to measure 7?.. 

The process forms a background when the tt" and 
e"*" interact in the target without leaving a detectable 
trace. Positron interactions were well-modelled in our 
EGS4-based simulation and we used the tt" energy 
deposition spectrum in scintillator measured previously 
in E787 [17| to model 7r~ absorption. We assessed the 
systematic uncertainty in the Kf,4, background by vary- 
ing the threshold of cuts on the energy deposited in the 
target fibers at the time of the pion. The kinematics cuts 
defining the signal region were 140 < P^r < 199 MeV/c, 
m < < 100.5 MeV and 12 < i?^ < 28 cm. We de- 
fined a sub-region 165 < P^ < 197 MeV/c, 72 < K < 
100 MeV and 17 < i?^ < 28 cm where the lower and up- 
per limits were chosen to suppress the K^i background 
that peaks near 160 MeV/c and the tail of the Kt^2 peak, 
respectively. 

The rejection of the ii'7r27 background was calculated 
using a combination of simulated if^2 and ^^,727 events 
and Kt^2 data events. The additional photon veto rejec- 
tion due to the radiative photon was calculated from the 
photon distribution in simulated events and the rejection 
power of single photons as a function of angle and energy 
evaluated with if^2 data [11]. 

Measurements of K^^ — > tt+tt" decay from the 
charge-exchange reaction were used as input to simulate 
charge-exchange events [Hj. The requirement on the de- 
layed coincidence between the reconstructed kaon and 
pion candidates provided suppression of charge- exchange 
background as the emitted tt^" was required to originate 
within the fiducial region of the target. The systematic 
uncertainty was assessed with the same methodology as 
the ife4 background. 

The muon and beam backgrounds were estimated en- 
tirely from data and were very small. As previous anal- 
yses had shown the muon background to be small 0] , 
the transient-digitizer-based cuts on 7r+ /i+ ^ e"*" 
identification were loosened to gain about 10% in ac- 
ceptance. The total acceptance of the signal region was 
(1.37 ±0.14) X IQ-^. 

In order to explore and verify the reliability of the back- 
ground estimates, we examined three distinct data re- 
gions just outside the signal region by loosening the pho- 
ton veto {PVn) or target pulse-shape (CCDn) cut. Each 
of the two regions, PVi and CCDi, were immediately 



Region 


Ne 


No 


P{No;Ne) 


Combined 


CCD^ 
PVi 
PV2 


0.79iiJ;* 
9.09ll;g 
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34 


0.45 [0.29,0.62] 
0.02 [0.01,0.05] 
0.61 [0.05,0.98] 


NA 

0.05 [0.02,0.14] 
0.14 [0.01,0.40] 



TABLE II: Comparison of the expected (A'e) and observed 
(A^o) number of background events in three regions CCDi, 
PV\ and PV2 outside the signal region. The central value of 
Ne is given along with the combined statistical and systematic 
uncertainties. P{No; Ne) is the probability of observing A^o 
events or fewer when A'e events are expected. The rightmost 
column "Combined" gives the probability of the combined 
observation in that region and the region(s) of the preceding 
row(s). The numbers in square brackets are the probabilities 
re-evaluated at the upper and lower bounds of the uncertainty 
on Ae M. 



adjacent to the signal region while a third region PV2, 
adjacent to PVi, was defined by further loosening of the 
photon veto cut. The number of expected and observed 
events and the probability of the observation are given 
in Table |TT1 The 5% probability for the regions nearest 
the signal region may have indicated that the background 
was over-estimated. Given the inability to cleanly isolate 
each background component by cut inversion, some con- 
tamination (i.e. events due to backgrounds from other 
sources) is possible and would generally inflate the back- 
ground estimates. Re-evaluation of the probabilities at 
the lower limit of the systematic uncertainties p^] gave 
14% for the two closest regions and demonstrated that 
the assigned systematic uncertainties were reasonable. 

After completion of the background studies, the signal 
region was examined and three candidates were found. 
The energy vs. range for these observed candidates is 
shown in Figure [T] along with the results of previous 
E787 d, 0] and E949 0,11] analyses. From these three 
new events alone, B{K~^ tt+i/P) = (7.89^ 
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X 10-^^ 
assum- 



was calculated using the likelihood method 
ing the SM spectrum and taking into account the un- 
certainties in the background and acceptance measure- 
ments [19j . When combined with the results of previous 
E787 and E949 analyses, we found B{K+ tt+vv) = 
{l.Tit\ll) X 10"^°- The signal-to-background (S/B) ra- 
tios for the three events are 0.20, 0.42 and 0.47 [13], 
which can be compared with the S/B = 0.20 for the pre- 
vious pnn2 candidate @ and with the S/B = 59, 8.2 and 
1.1 for the pnnl events 0] assuming B{K^ ^ ■k^vv) — 
1.73 X 10^^°. In this analysis, a candidate in the best 
(worst) cell would have had S/B=0.84 (0.20). The prob- 
ability that the three new events were due to background 
only, given the estimated background in each cell, is 
0.037. The probability that all seven K'^ — > -k^vv events 
were due to background is 0.001. In summary, these ob- 
servations imply a -k^vv branching ratio consis- 
tent with SM expectations. 
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FIG. 1: Kinetic energy vs. range of all events passing all other 
cuts. The squares represent the events selected by this anal- 
ysis. The circles and upward-pointing triangles represent the 
events selected by the E787 and E949 pnnl analyses, respec- 
tively. The downward-pointing triangles represent the events 
selected by the E787 pnn2 analyses. The solid (dashed) lines 
represent the limits of the pnnl and pnn2 signal regions for 
the E949 (E787) analyses. Despite the smaller signal region 
in Et^ vs. Rtt, the pnnl analyses were 4.2 times more sensi- 
tive than the pnn2 analyses. The points near i?^ = 108 MeV 
were /(',r2 decays that survived the photon veto cuts and were 
predominantly from the pnnl analyses due to the higher sen- 
sitivity and the less stringent photon veto cuts. The light 
gray points are simulated if^ -k^vv events that would be 
accepted by our trigger. 
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